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1 Introduction 

Reduction by symmetry of Lagrangian field theories has aided the implementation of many 
diverse mathematical models from geometric mechanics. Two main approaches have been 
developed. One approach, investigated in Gotay et al. [2004], employs multisymplectic ge- 
ometry to extend the symplectic formulation of classical Lagrangian systems. The second 
approach, studied in Castrillon-Lopez et al. [2000] and Castrillon-Lopez & Ratiu [2003], re- 
duces the variational principal itself without reference to the Hamiltonian side and is referred 
to as covariant Lagrangian reduction. 

A comparison of the covariant Lagrangian reduction approach, Castrillon-Lopez et al. 
[2000]; Castrillon-Lopez & Ratiu [2003], with the corresponding classical Lagrangian reduc- 
tion method. Holm et al. [1998]; Cendra et al. [1998, 2001], shows that a paradigm permitting 
both reductions is currently lacking. Further, examples such as the isoperimetric problem 
and metamorphosis image dynamics require such a paradigm. And further still, the desired 
capability for covariant reformulations of classical problems, such as in Marsden & ShkoUer 
[1999], call for such a paradigm. 

This paper achieves a full generalization of the classical theory while preserving the flavor 
of the current covariant theory. Applications to the Skyrme model and the molecular strand 
illustrate the ideas of CastrilkSu-Lopcz et aL [2000] and CastriUon-Lopcz & Ratiu [200-)] in 
§1.1 and §1.2, respectively. A general discussion of classical Lagrangian reduction appears in 
§1.3. These discussions illustrate the need for the development of the more general theory. 

1.1 Principal bundle reduction: the Skyrme model 

The first in the series of papers on covariant Lagrangian reduction, Castrillon-Lopcz et al. 
[2000], dealt with the extension of classical Euler-Poincare reduction of variational principles 
to the field theoretic context. There, a field theory was formulated on a principal bundle 
and was reduced by the structure group. These results may be illustrated by the Skyrme 
model for pion interaction, which was first developed in Skyrme [1961] and whose more recent 
developments were reviewed from a Physics-based standpoint in Gisigcr & Paranjapc [1998]. 

Both the original formulation of the classical Skyrme model and its recent advances have 
been described, for example in Gisiger & Paranjape [1998], in terms of local coordinates. 
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The use of local coordinates, while necessary for numerical implementation, tends to obscure 
the geometric content of the equations. Therefore this paper approaches the theory, on the 
whole, from a coordinate-free viewpoint. However, to aid communication and compatibility 
with the references, some of the examples are addressed in local coordinates too. 

The Skyrme model. The class of model given in Gisiger & Paranjape [1998] may be 
outlined as follows: 

Consider a unitary field u : ^ SU (n) over the three-sphere 5*^ for n either 2 or 3. 
The three-sphere is interpreted as a one-point compactification of Euclidean space M^. In 
local coordinates, the massless Skyrme Lagrangian reads 

Lsky {u,dju) = -^tr {u^ {dju)u\dju)) + T^tr (^[u\dju),u\dku)Y^ (1.1) 

where := v/^ is the adjoint of u. The constants and are potentially calculable 
from QCD but, in practice, are fitted to experimental data. The local representation of the 
Euler-Lagrange equations for Lsky are given by 

d, (^Y^d^n)u^ + ^ [{dku)u\ [{d^u)u\ {d'u)u^]]^ = 0. (1.2) 

Baryons are identified with topological soliton solutions of equation (1.2) with u : 
SU{n). Note that the Lagrangian Lsky is S'[/(ri)-invariant under the transformation 

{u,dju) {vu,v{dju)) , for all v G SU{n). 

Therefore, a reduction by symmetry may be effected. In order to bring out the geometry of 
the system a reformulation of the problem is required. 

Geometric formulation. Let vr := tts'-^.p : -P — > 5*^ be a principal SU{n) bundle over the 
three-sphere S^. A section of vr is a smooth map a : ^ P, such that 

n o a = ids3, (1.3) 

where id^a is the identity map on S^. The space of sections of vr is denoted T (vr). Recall 
that a principal bundle admits a section if and only if it is trivial. Therefore, in general only 
local sections a : U C ^ P defined on an open subset U C may be considered. 

In a local trivialization U C 5''^, a section a reads a{x) = {x,u{x)), where u : U G ^ 
SU{n). Thus, the space of sections T (tt) corresponds to the space of unitary fields. 

Recognize that {u, dju) is a local representation of the tangent map Ta : TS^ — *• TP. 
The jet bundle, J^P, provides the natural space to consider such objects. This affine bundle 
over P may be defined fiberwise by 

j;P = { 7p G {T.S^ TpP) I Tpvr o 7^ = idr^s^ } , 

with projection vrp^ip : J^P P given by 7rp jip(7p) = p. The jet bundle serves field 
theories as the tangent bundle serves classical Lagrangian systems. 
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For the most part, J^P is considered as a fiber bundle over with projection 

TTgs ji p := n o 71 p ji p : J^P — ^ . 

Indeed, the tangent map of a section a G r('7r), interpreted as a map x ^ TxCr, offers a 
section of ns^^jip since for T^cr G L (T^S^ ,T„(x)P) equation (1.3) yields 

T^(^)7r o T^a = (vr o a) = TJdss = idx^s^- 

The geometry introduced here is succinctly visualized and organized by commutative 
diagrams. The following commutative diagram exhibits the geometry of the jet bundle: 

jip ^If^ . 53 




P 



Arrows between spaces indicate maps from the space at the tail of the arrow to the space at 
its head. Sometimes arrows are adorned with the name of the maps they represent. Different 
paths through the diagram are equivalent in terms of composition of the associated maps; 
therefore, this diagram also communicates the relation 

7rs3,JiP = O TTpjip. 

The reduced bundle. Having identified the geometry of the Classical Skyrme Model, one 
may proceed by thinking about reduction by left SU{n) symmetry in the style of C'astrillou- 
Lopez ot al. [2000]. The quantities {dju)u^ may be understood as the local representation 
of a principal connection form A G Q^{P,su{n)) which has been pulled back by the unitary 
field where su{n) denotes the Lie algebra of SU{n). The connection form A G Q^{P,u{n)) 
on P provides the required geometric tool to effect the reduction since it provides a vector 
bundle isomorphism 

J^P/SU{n) L (T5^ AdP) , (1.4) 

where AdP denotes the adjoint bundle associated to the principal bundle P defined as the 
quotient space 

AdP := (P X su{n)) /SU{n), 
relative to the following diagonal action of m G SU{n): 

(p, G P X su{n) 1-^ {up,AduO £ -P X 3u{n). 
Denoting the equivalence class of (p, ^) G P x su{n) by 

b,a.u(n)eAdP, 

the bundle isomorphism J^P/ SU{n) —>■ L (TS*^, AdP) reads 
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The reduced Euler-Lagrange equations. The Skyrme model is now written in the 
same form as the resuh of . ... .rillou-Lopez et al. [2000] which states that the Euler-Lagrange 
equations on J^P are equivalent to the covariant Euler-Poincare equations on L {TS^, AdP), 
which read 

div^|i-ad;|^ = 0. (1.5) 
oa OCT 

Here div"^ denotes the divergence associated to the covariant exterior derivative on AdP 
associated with the principal connection A and ad* is the dual of the adjoint operator on 
su{n). For the classical Skyrme model, the reduced Lagrangian associated to (1-1) can be 
written as 

where || ■ || is the norm associated with a Riemannian metric on and an Ad-invariant inner 
product on g. The unreduced Lagrangian density is 

and is clearly S'f/(n)-invariant. The classical Skyrme model equations then become 



div-^n - ad;n = n = + A^ad; (a a c^)^ (1.6) 

2 8e^ 



where b denotes the flat map b : L{TS^, AdP) —>■ L{T*S'^, Ad*P) induced by the Riemannian 
metric on 5*^ and the Ad-invariant inner product on g. The local representation of equation 
(1.6) is equation (1.2). For details on related dynamical systems to the classical Skyrme 
model see Holm [2008] . The link between the covariant and dynamical reductions associated 
to the equation (1.5) is established in Gay-Balmaz & Ratiu [2009]. 



1.2 Subgroup reduction: the molecular strand 

The Skyrme model illustrates reduction of a principal bundle by its structure group as de- 
scribed in Castrillon-Lopez et al. [2000]. Correspondingly, the molecular strand demonstrates 
reduction of principal bundles by a subgroup of the structure group, which was the subject 
of C'astrillon-Lopcz & Ratiu [2003]. A molecular strand may be modeled as a flexible, elastic 
filament moving in 

with rigid charge conformations undergoing rigid rotations mounted 
along the filament's length, as shown in Figure 1.1. A full treatment of the molecular strand 
was undertaken in Ellis et al. [2008]. 



The geometry of the molecular strand. The parameter space for the molecular strand 
is X = / X M where / is an interval of fixed length and M represents time. The strand may 
be described by maps 

(A, r):X^ SE{3) ^ S0{3) x x= {s, t) ^ (A(x), r(x)) . 

Here r : X — > describes the position of a point on the filament at a given time and 
A : X — >■ 50(3) describes the rigid charge conformations along the filament at a given time. 
These maps correspond to sections a of the principal bundle 

7rx,P : P := X X 5E(3) ^ X, 7rx,p(x, A, r) := x, 
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Figure 1.1: Rigid conformations of charges are distributed along a curve, 
by the relation 

cr(x) = (x, r{x), A(x)) . 

The Lagrangian L : J^P ^ M of the molecular strand is taken to be left S'0(3)-invariant 
as in Ellis et al. [2008]. Contrary to the case of the Skyrme model, the symmetry group 
of the theory does not coincide with the structure group SE{3) of the principal bundle P, 
but is a is a subgroup of SE{3). Thus, departing from the Skyrme model, P is a principal 
5'0(3)-bundle with the projection 

TTs.p : P ^ S := X X vrg^p (x, r. A) = {x, A" V) . 

Now there are two bundle structures on P given by 7rx,p and tc-^^p- These two bundle 
structures induce a third, this time on S = X x M^, given by 

TTx.s : S X, ttx,e{x, p) = x. 
The geometry is described by the commutative diagram. 

X X SE{3) X 




X X M3 



The reduced bundle. Since the symmetry group is a subgroup of the structure group, 
there is an additional part to the bundle isomorphism (1.4). More precisely, a principal 
connection A on 7rs,p furnishes J^P/S0{3) with the vector bundle isomorphism 

J^P/S0{3) J^E ©E L (tt^ s^^, AdP) 
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given by 

[Tcr]soi3) ^ {t (7rs,P o a) , {a, a*Aj^^^.^^^ . (1.7) 

This isomorphism will be studied in detail in §2.2. Since, for the molecular strand, vrs,p is a 
trivial bundle, the adjoint bundle Ad P is also trivial and can be identified with E xso(3) S 
(projection on the first factor) via the isomorphism [(a;. A, r), ^J^^^g^ i— >• (x, A~^r, AdA-i ^). 
Using the Maurer-Cartan connection 

the isomorphism (1.7) reads 

[(x, r, A, id, dr, dA)]gQ(^^'^ ^-^ (x, p, dp, A^^dA) = (x, p, pdt + p'ds, tudt + fids) 

where p = A~V, u = A^^A, Q = A^^A', (') denotes differentiation with respect to t, and 
(') denotes differentiation with respect to s. 



The reduced Euler-Lagrange equations. The main result from Castrillon-Lopcz & 
Ratiu [2003] states that when nx,p : P ^ X is a principal bundle and G is a subgroup of the 
structure group, then the Euler-Lagrange equations on P for a left G-invariant Lagrangian 
L are equivalent to the Lagrange-Poincare equations for the reduced Lagrangian / on the 
reduced bundle J^P/G ^ J^HQj^L [ti^^^TX, kdP) . The La grange-Poincare equations read 

_4 51 61 
div — - ad^— = 0, 

da da , . 

diviL _ - 
5dp 6p 

The exact definitions of the differential operators and the functional derivatives appearing 
in these equations are studied in §2 and §3. The right hand side of the first equation in (1.8) 
usually has an extra term associated with the curvature of the principal connection A. This 
term is described in more detail below, but for the present case A is fiat, so the curvature 
term vanishes. In local coordinates equations (1.8) read 

ouj d\l 
6p ^5p' 6p 

These equations need to be augmented with an integrability condition to allow reconstruc- 
tion. This integrability/reconstruction condition is related to the curvature of A. For the 
molecular strand the required reconstruction condition is 

(9jfi — dgUJ — u X Q = 0. 

These equations recover the results derived in Ellis et al. [2008] where the Lie algebra, so (3) 
was identified with M.^ and therefore the adjoint actions became cross-products. More details 
about the reconstruction condition can be found in 53.3. 
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1.3 Fiber bundle reduction 



In §1.1 we introduced the ideas of Castrillon-Lopez et al. [2000] in the context of the classical 
Skyrme model. Reduction of such 'pure gauge' theories requires the introduction of certain 
geometric tools such as the adjoint bundle and jet bundles. In §1.2 we used the example of 
the molecular strand to review how principal bundle reduction can be extended to include 
reduction on principal bundles by a subgroup of the structure group as in Castrillon-Lopcz 
& Ratiu [2003]. The purpose of this paper is to extend these ideas farther to include fiber 
bundle reduction. 



Lagrangian reduction in classical mechanics. Consider classical Lagrangian reduction 
from the perspective of reduction of variational principles. A variational principle is formu- 
lated on a principal bundle 7i : Q ^ Q/G and a principal connection A is introduced on Q. 
The connection yields a bundle isomorphism 

(TQ)/G^TiQ/G) ©Q/cAdg 

given by 



[v,]^ (T7riv,),lq,Aiv,)l 

Thus, a curve q(t) E Q induces the two curves 

p{t):=niqit))eQ/G and = |g(t), ^ (g(t))]g G AdQ. 

Classical Lagrangian reduction states that the Euler-Lagrange equations on Q with a G 
invariant Lagrangian L are equivalent to the Lagrange- Poincare equations on TQ/G = 
T{Q/G) ©Q/c AdQ with reduced Lagrangian /. The Lagrange-Poincare equations read 




where B is the reduced curvature form associated to the connection A and D/Dt denotes 
suitable covariant derivative. 



Extension to field theories. When this classical approach is generalized to field theories 
the natural structure to consider is the trivial fiber bundle 

7rx,p : P := X X Q ^ X, 'itx,p{x, q) := x. 

Now r(7rx,p), the space of sections of vrx,p, generalizes the space of curves in Q and the 
principal bundle structure vr on Q gives a natural principal bundle structure 

TTs.P :P:=Xxg^S:=Xx {Q/G) , 7rE,p(x, q) := {x, 7r(g)). 
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More generally, one would consider the following commutative diagram 



P 




where 7rx,p : P ^ X is any fiber bundle and vr^^p : P ^ S is a principal bundle, whose group 
action preserves the fibers of 7rx,p- This situation arises, for example, in the isoperimetric 
problem and in image dynamics, as is outlined in §5. 

Goals of the paper. Following the preceding discussion, a specification of properties 
emerges. One would like to develop a framework for reduction that: 

1. Captures the natural generalization of classical Lagrange- Poincare reduction to field 
theories. 

2. Reduces to classical Lagrange-Poincare reduction as a particular case. 

3. Reduces to the work of Castrillon-Lopez et al. [2000] and Castrillon-Lopez & Ratiu 
[2003] in particular cases. 

In §2 some geometric tools that are necessary for performing reduction are introduced. 
The relationship between two bundle structures on the same manifold P is also studied. In 
§3, the Lagrange-Poincare field equations are developed and a method for reconstruction 
is given. The Kelvin-Noether theorem associated to the Lagrange-Poincare field equations 
is presented in §4. Finally, in §5, the reduction tools developed earlier are applied to the 
isoperimetric problem and image dynamics. Throughout the paper, constant reference to 
four particular cases is made in order to illustrate the theory: the unreduced case, principal 
bundle reduction, subgroup reduction, and classical reduction. 

2 Geometric constructions 

There are two main geometric constructions of interest. The first is the interaction of two 
bundle structures, 7ix,p and tts p on P. The second is the reduction of the jet bundle J^P 
by the structure group G. 

2.1 Geometric setting 

Consider a locally trivial fiber bundle 7rx,p : P —>■ X. A section of ttx,p is a smooth map 
a : X P such that t^x,p o a = idx- We denote by T{nx,p) the set of all smooth sections 
of 7rx,p. 
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Remark. It is necessary to introduce many fiber bundle projections during the develop- 
ment of the theory. The notation indicates the source and target space, e.g. vrx,p : P — >■ X, 
where the first subscript denotes the base space and the second the total space. The order 
of the subscripts allows one to write, for example. 

The first jet bundle of 'itx,p is the affine bundle vrp jip : J^P P whose fiber at p is the 
afiine space 

JpP = {7p e L{T^X, TpP) I Tpnx,p oip = idr^x, x = nx,p{p)} , 

where L{TxX,TpP) denotes the space of linear maps 7p : T^X — > TpP. The first jet bundle 
is the natural generalization of the tangent bundle to the field theoretic context. Therefore 
J^P plays the role of our unreduced state space in applications. The manifold J^P may 
also be regarded as a locally trivial fiber bundle over X, that is, t^x,j^p '■ J^P ~^ X with 
^x,jip(7p) '■= '^x,p{p)- Given a G r(7rx,p), the first jet extension of a, defined by jV(a;) : = 
Ti,.cT, is a section of tcxj^p- 

Suppose there is a free and proper left action $ of a Lie group G on P such that 

vrx,p °^g = vTx.p, for all g e G. (2.1) 

Equation (2.1) is equivalent to the assumption that the action of G preserves the fibers of 
7rx,p. Since the action is free and proper, there exists a principal bundle n^ p : P S, 
where S := P/G. Here S is the equivalent of shape space in applications. Since, by (2.1), 
the projection itx,p is G- invariant, it induces a surjective submersion 7rx,E : S — > X via the 
relation 

7rx,E o 7rs,p = vrx,p. (2.2) 

It is easily verify that if 7rx,p is proper then 7rx,s is also proper. 

More generally if 7rx,p is a locally trivial fiber bundle then vr^.s is also a locally trivial 
fiber bundle. To see this, take a fiber bundle chart ip : 7r^^p{U) ^ U x F, where U is an open 
subset of X and the manifold F is the model of the fiber. By definition, pioif) = 7Tx,p, where 
Pi : U X F ^ U is the projection onto the first factor. Property (2.1) implies that 7Cx^p{U) 
is a G-invariant subset. Thus, the diffeomorphism ip bestows a well-defined G-action on 
U X F which turns out to be free, proper, and acts only on the component F. The model 
fiber F thereby attains a principal bundle structure F — > F/G induced by (and depending 
upon) the chart ip. Since ip is an equivariant diffeomorphism, it drops to a diffeomorphism 
ijj : 7r^p([/)/G ^ U X F/G. Also, since ■K^p{U) / G = 7r^^j.{U) and pio ijj = 7rx,s, the map 
is a fiber bundle chart of ttx.e : S — > X. For principal bundles, one needs to work with 
local sections, since a principal bundle does not have global sections unless it is trivial. 

Remark. There are now two different bundles tix,p '. P —>■ X and 7rs,p : P —>■ with 
the same total space P. In general, the associated vertical distributions do not coincide, 
although (2.1) provides the inclusion 

ker (Tp7rs,p) C ker (Tp7ix,p) ■ 
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Thus, it is possible to associate two different jet bundles to P. Throughout this paper only 
the first jet bundle of interest is 7rx,p : P ^ X and hence there is no ambiguity in the 
notation J^P. 

Lagrangian field theories are described by a Lagrangian density C : J^P A"+^X 
defined on the first jet bundle. Here K^^^X denotes the bundle of (n + l)-forms on X, where 
n + 1 = dimX. In this context the G-action on P, lifted to J^P, should be interpreted as 
a symmetry of the Lagrangian density. The associated reduction process, described in the 
next section, is called the covariant Lagrange-P oincare reduction. 

Particular cases. Various previous theories may be identified as particular cases of the 
geometric setting developed in this paper. These examples will be referred to throughout 
the paper and serve to illustrate the ideas introduced in more familiar contexts while demon- 
strating how the objective of capturing previous theories in this new context is fulfilled. 

i If S = P, that is, G = {!}, there are no symmetries. The principal bundle structure 
disappears and the geometric setting for covariant Lagrangian field theory, referred to 
as the unreduced case, emerges. The commutative diagram that describes this case is 



P ^X 




E = P 



where vr^ p is a fiber bundle. 

ii Assume that the configuration space 7rx,p : P ^ X is itself a principal G-bundle 
and G is also the symmetry group. Then S = X and vrx,E is the identity map. 
This recovers the geometric setting for covariant Euler-P oincare reduction, or principal 
bundle reduction in Castrillon-Lopez et al. [2000] which is used to study, for example, 
the Skyrme model. The commutative diagram that describes this case is 



P ^X 




S = X 



where nx^p is a principal bundle. 

iii If Tix^p : P ^ X is a principal bundle whose structure group H contains the group of 
symmetries G as a subgroup one recovers the formulation in Castrillon-Lopez & Ratiu 
[2003]. This is the geometric setting for the formulation of the molecular strand from 
Ellis et al. [2008] and is referred to as subgroup covariant Lagrange-P oincare reduction 
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or simply subgroup reduction. The commutative diagram that describes this case is 



P 




E 



where 7rx,p and vr^^p arc respectively G and i7-principal bundles. 

iv If X = M, S = MxM, and P = RxQ, where tt^-^q : Q ^ M is a G-principal bundle, the 
geometric setting for Lagrangian reduction in classical mechanics, known as classical 
reduction, becomes apparent. Here Q plays the role of the configuration space. There 
are two well-known particular cases: li Q — G (thus M — {m}) the geometric setting 
for Euler-Poincare reduction surfaces (this is also a particular case of ii), where the 
configuration space coincides with the group of symmetries; if G = {1} (thus Q = M) 
there are no symmetries and we reacquire the geometric setting for unreduced classical 
Lagrangian mechanics (this is also a particular case of iii) . The commutative diagram 
that describes this case is 

P^RxQ "If^i^^ -X = K 




E = M X M 



where ttm.q is a principal bundle. 



Adjoint bundle. The adjoint bundle associated with the principal bundle 7rs,p : P — > E 
is a vector bundle 7rs,AdP '■ AdP — > E. The total space AdP is the quotient space P Xq Q 
relative to the following left action of G on P x 

(p,0^($p(p),Ad,0- 

Elements in the adjoint bundle are equivalence classes [p, ^J^ and the projection is described 

by 

7rs,AdP (b>Clfl) = 7rs,p(p). 

The adjoint bundle is, in fact, a Lie algebra bundle. That is, each fiber (AdP)s, s e E, has 
a natural Lie bracket 

where 7rs,p(p) — s, ^,r) & Q and these Lie brackets depend smoothly on the base variable 

s e E. 
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This Lie algebra bundle structure enables the introduction of a wedge product. For 
1-forms this wedge product A : ^^(S, AdP) x ^]l(E, AdP) ^ n^{J:,AdP) is defined by 

{a A p){us, Vs) ■■= [a{us),Pivs)]s - [a(f,), (2.3) 

where s G S and Us,Vs G T5E. 

The different equivalence classes are interpreted as different representations of the dynam- 
ics. For given s G S and p G vr^ p (s) one can define a p-dependent Lie algebra isomorphism 
Xp : (AdP), ^ g by 

Xp(lq,vL)=t (2.4) 



where ^ G 3 is such that (p, = [q, r]}^. 

The choice of p G tt^ p (s) determines the representation of the dynamics. Thus by 
altering p we use Xp to give the dynamics in the convective or the spatial representation. 

A connection A on the principal bundle vr^.p : P — S is a one-form A G fi^(P, g) such 
that 

<i>lA = AdgoA and A{^p{p))=^, where ^p{p) := - $exp(to(p) 

at j^g 

is the infinitesimal generator associated to the Lie algebra element ^ G g. The horizontal 
distribution associated to A is the subbundle HP C TP defined by 



HpP:=keT (Aip)). 

The horizontal distribution is complementary to the vertical distribution VpP := ker (TpTr^p) 
and consequently TP decomposes as TpP = HpP © VpP. The connection A defines the 
horizontal lift operator Hor^ : T,!! HpP according to 

B.OTp{vs) G HpP and T^ttep o Hor^ = idr.s 

where s G S, p G vr^^ p (s) and Vs G T^E. 

The connection A also induces a covariant derivative on the adjoint bundle 

V-^ : r(7rs,Adp) r (7rs,L(TS,AdP)) (2.5) 

given by 

a{c{t)) = lp{s),davs) - [A{Tp{v,)),as)]i, (2.6) 



t=0 

where ^ : S — > g and p : T, ^ P are such that cr(s) = |p(s),^(s)]g and c(t) is a curve in S 
such that c(0) = (see Cendra et al. [2001], Lemma 2.3.4). The covariant derivative V"^ 
also has an interpretation bilinear map 

V-^ : X(S) X r (7rs,Adp) ^ F (vrs.Adp) , (X, a) ^ Via. 
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2.2 Reduced covariant configuration space 

The free and proper action $ : G x P ^ P induces a free and proper action $ 
J^P defined by 

<l>;(7p):=T$,o7„ ^,e4P. 
Note that this action preserves J^P since, by (2.1), 

Tnxp o $g(7p) = Tnxp o T^g 07^ = Tttxp o 7p = idr^x- 

Thus it is vahd to consider the quotient manifold J^P/G 3 [1p]g- 

Remark. Recall that J^P denotes the first jet bundle of P as a fiber bundle over X and 
not clS cl principal bundle over S. 

The connection A on the principal bundle ttsp : P ^ S introduces the smooth map 
which is defined by 

(3^ : jip/G ^ jiS©s^ (7ri_2TX,AdP) , /3^([7Jg) := (T7rs,p07p, b, ^(7p(-))lg), (2-8) 
where 

^s,L(.i_^TX,Adp) : L{t:\ .^TX,K<\P) S 

is the vector bundle whose fiber at s G S is L(T7r^ Ad Ps). 

The map is an diffeomorphism, the inverse being given by 

(3-/ ■ J^T.®^L{7T\^^TX,KdP) ^ J^P/G, f3~/{5M = [Hor^ o 5, + Cp o /,] ^ , 

where p G P is such that vr2,p(p) = s and 

Cp:{MP)s^VpP 

is defined by Cp([Q')'7]g) = ^p(p) where G g is such that If?, "^^Ig = |p, ^]g- Note that 
Hor^ o5s + Cp°h £ J^P and that its equivalence class does not depend on which p is chosen. 
Also, note that the diffeomorphism /3_4 endows the manifold J^P/G with the structure of an 
affine bundle over S. 

The isomorphism is interpreted as follows in the four particular cases: 

i Here G = {1} thus the principal bundle structure disappears. The bundle isomorphism 
is the identity on J^P. 

ii Here S = X, thus j3_A is a bundle map over X and we have 

; jip/G ^ L(TX,AdP), /3^([7p]) = b,^(7p(-))lg 

and we recover the isomorphism used in the Euler-Poincare reduction, see formula (2.5) 
in Castrillon- Lopez et al. [2000]. 

iii The isomorphism used in the subgroup Lagrange-Poincare reduction reemerges, see 
Proposition 3 in Castrillon-Lopcz & Ratiu [2003]. 



: Gx Jip ^ 
(2.7) 
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iv Since X = M and P = M x g, the jet bundle J^P may be identified with M x TQ. 
Thus, J^P/G ~ M X (TQ/G). Similarly, J^S may be identified with M x TM and 
AdP with M X AdQ. The bundle J^E ©s L {ti*x j.TX, AdP) can thus be identified 
with M X (TM (Bm AdQ). A connection 7 on Q naturally induces a connection A on 
P and the bundle map (3a reads : M x (TQ/G) ^ M x (TM ©m AdQ), 

/5.4(t, K]g) = {t,TnMM, [g,7(^g)lg) • (2.9) 

Therefore the usual connection dependent isomorphism TQ/G ~ TM ©^ AdQ used 
in classical Lagrangian reduction is recovered, as in Ccndra ct al. [2001]. 



3 Lagrange- Poincare field equations 

Consider a G-invariant Lagrangian density C : J^P A"+^X. For simplicity, suppose that 
X is orientable and fix a volume form /i on X. The Lagrangian density may thereby be 
expressed as £ = Lfi, where L : J^P M. 

Let U C X be an open subset whose closure U is compact. Recall that a section 
(T : [7 — > P of 7i'x,p is, by definition, smooth if for every point x & U there is an open 
neighborhood of x and a smooth section '■ ^ P extending a. A critical section of 
the variational problem defined by C is defined as a smooth local section a -.tj P oi 7ix,p 
that satisfies 



de 



0, 



e=0JU 

for all smooth variations a^, -.IJ ^ P such that = a and (Te\du 



5a{x) :- 



d 
de 



(7\du- Since 

ae{x) G K(x)-P := ker(T^(^)7rx,p) and (5cT|ac/ = 0, 



£=0 



one may assume without loss of generality that = 0^ o a, where 0^ is the fiow of a vertical 
(with respect to the bundle structure vr^p) vector field V G X^(P) such that V{a{x)) = 
for all X G dU. The smooth Tietze extension theorem facilitates V^s definition over the 
whole manifold P, but values of V outside cr{U) will not play any role in any subsequent 
consideration. Note that Sa{x) = V{a{x)) for all x E Lf. Consequently, 



(5jV(x) := 



d 

de 



jVe(x) = jV(jV(x)) G Vji^(^) J^P = ker(T,i^(^)7rx,jip), 



£=0 



where V G 3e^(P) ^ j^V G X^(J^P) is the 1-jet lift of vector fields. Thus, 0" is a critical 
section of the variational problem defined by C if 







d 

de 



=0 Ju 



5L 



x),jV(jV(x)))/i 



where {6L/6a) (x) G ^*o-(i.)<^^-P is the differential along j^a. That is, 

/6L 



— (x),Z(jV(x)) > = dL(jV(x))(Z(jV(x))) 
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for arbitrary vector fields Z e X^(J^P) that are vertical with respect to 'Kx,j^p- Denoting 
by ELiV) the bundle morphism £C{L) : J^P V*P defined by the condition 

^(6:£(L)(jV(x)),r(a(x))>/i = ^^^(x),jV(jV(x))^^, for all V G X^(P), 

the covariant Euler-Lagrange equations can be written intrinsically as 

SC{L) = 0. 

Here SC is represented locally by 



£C{L) 



1 d f dL , 



where C = L{x'^ , y"' , v°')d"'~^^x . Thus, in coordinates, the covariant Euler-Lagrange equations 
take the standard form 

'^A)-#.fl^(A)l=o. (3.1) 



dy"' dx^ \dvi 

These equations may be written globally by using a connection on the affine bundle iTpjip : 
J^P — > P; this point of view will be used at the reduced level. 

By G-invariance, L induces the reduced Lagrangian I : J^P/G M. Fixing a connection 
A on the principal bundle vr^^p : P ^ S brings in the bundle isomorphism thereby 
permitting the definition of the reduced Lagrangian / on J^E ©E L {■K^ j.TX, Ad P) . 

A section a G V{-kx,p) of the configuration bundle induces a section 

p ■= T^T.,p OCT e r(7rx,s) 
by (2.2). The reduced section is defined as 

a := la,a*Al E T (vr^,^(.^^^^,AdP)) • (3-2) 

Thus, 

(jV,^) = Pa ([jV]g) : X ^ J^S ©s ^ (vr^^^TX, AdP) . 
The two components are not independent since p can be obtained from a; explicitly. 

Note that {j^p,a) is a section of the bundle J^S ©^ L (vr^^^ ^^TX, Ad P) viewed as a fiber 
bundle over X, and not as an affine bundled over S. These definitions and the G-invariance 
of C (and hence of L) yield 

£(jV)=L(jV)^ = /(jV,fT)/x (3.3) 

for any a G r(7rx,p). 

The previous considerations hold without changes when a is a local section a : U ^ P. 

The fact that J^S ©s L {t^\^^TX,MP) is a locally trivial fiber bundle over X follows 
from the following observation: G acts on the locally trivial fiber bundle 'Kx,j^p '■ J^P X 
by a free and proper action such that vr^.jip o = 7rx,jip. Therefore, by the argument 
used in §2.1, J^P/G — >■ X is a locally fiber bundle. Thus, the isomorphism ensures that 
J^E ©2 L (^71^ X , Ad P^ is a locally trivial fiber bundle over X. 



16 



3.1 Reduced variations 



Using the bundle isomorphism /3_4, the variation of the action defined by L gives 

e=0JU 

where ?7 C X is an open subset and cTe : U —>■ P is a. smooth variation of the smooth section 
a:U ^ P. 

A covariant derivative on the locally trivial fiber bundle 7!'x,AdP '■= 7rx,s°7i"E,AdP : Ad P — *• 
X is required to compute the reduced variations. Recall the following general construction. 



Ju 



General constructions. Let r S a vector bundle endowed with a covariant deriva- 

tive V. Recall that V induces a covariant exterior derivative : f2'^(E,£') — > 
whose formula is a direct adaptation of the standard Cartan formula for fc-forms on a man- 
ifold, by replacing all directional derivatives by covariant derivatives relative to V. In par- 
ticular, for one-forms 

d^aiU, V) = VuHV)) - Vv(a(f/)) - V^]), (3.4) 

where a G E) and U,V e 

Let X be an arbitrary manifold and f : X ^ E a. smooth function. Define the V- 
derivative of / by 



/(c(()) £ (3.5) 

t=0 

V 



where c{t) G X is a curve such that c(0) = and D /Dt is the usual covariant time deriva- 
tive associated to V of the curve t i— >• /(c(t)) in E. Note that V/(x) e L (T^X, Er{f{x))) , 
and when U G X(X), the derivative V[// is a function on X taking values in E. 

The considerations below require an exterior covariant derivative of forms on X with 
values in E. To make sense of this, assume that there is a smooth map tt : S ^ X. Recall 
that r'= : L'^{tt*TX,E) ^ S is the vector bundle whose fiber at s G S is L'^{tt*TX, E)s = 
L^{TTr(^s)X, Eg), the fc-linear antisymmetric maps from T^(s)X to Eg. Define the i?-valued 
fc-forms on X by 

^ti^^E) ■= r {TTx,L'^{n*TX,E)) , 

where t^x,l'^{-k*tx,e) := vr o r'^ : L^{t:*TX, E) — > X. Note that this is not a vector bundle 
and thus ^^^(X, E) are not the usual vector bundle valued fc-forms on X. In fact, f2^(X, E) 
is not even a vector space. In spite of this, there is a derivation, analogous to the usual 
exterior covariant differentiation (3.4) on forms. While the definition of this operator holds 
for general elements in Vt^{X, E) and is again based on Cartan's classical formula, only the 
definition for one-forms is needed: 

d^aiU, V) = VuiaiV)) - Vv(a(f/)) - ai[U, V]), (3.6) 

where a G ^ll{X, E) and U,V e X(X). Note that a(f/), a{V) : X ^ E, hence (3.5) is vahd, 
also note that d'^a G Ql{X,E). 

Since Ll{n*TX, E) = E hy definition, = r : E ^ S and thus fi" (X, E) = T{7rx,E) C 
C°°{X,E), Trx,E := vr o r°. Therefore the operator d^ on fl^{X,E) coincides with V as 
defined in (3.5). 
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Covariant derivatives. Returning to the case at hand, the general construction special- 
izes to the covariant derivative V"^ on the vector bundle AdP — >• S. Thus if ^ G F (vrx.Adp) 
and G T^X, the previous definition, 3.5, gives the V'^-derivative of ^ by 



t=0 

where c(t) is a curve in X such that c(0) = u^. Writing ^(s) = C(2^)L yields 



Vie(a:) = dCK) - [A{TMu.)), ax)]}^ (3.7) 

(see Cendra et al. [2001], Lemma 2.3.4). Note that the V'^-derivative is a map 

V-^ : r (vTx.Adp) r (^1T x,LiTr*^^TX, Ad P)^ , 

not to be confused with (2.5), and it can be interpreted as a map 

V-^ : X(X) X r (7rx,Adp) ^ T (vr^.Adp) • 
Note also that V^^^, V;^^, and ^ project to the same section p G r(7rx,s), that is, 
^s,L(^^^rx,Adp) ° ^'^^ = 7rs,AdP o = 7rs,AdP o ^ = P- 
Next, in the present situation, the covariant exterior derivative 

:= d^^ : Ql^^^ (X, Ad P) ^ Ql^^^ (X, Ad P) (3.8) 
is attained from (3.6). That is, 

d''a^^.,v,) = m)) - m)) -a[u,v]), (3.9) 

where ^ G ri^^_^(X, AdP), u^,v^ G T^X, U,V e X(X) satisfying U{x) = u^, V{x) = v^. 

Variations. The wedge product (2.3) in ^^(S, AdP) extends to elements a, /3 G f2^^^(X, AdP) 
that project to the same element p G r(7rx,s) by 

(a A /3)(u^,fa;) := [a{ux),P{v^)]p(x) - [a^v^), P{ux)]p{x), (3.10) 

where Ux,Vx G T^X. 

Proposition 3.1 Let a : X ^ P be a smooth section of 7ix,p : P —>■ X . Let A be a 

connection on the principal bundle txy,,p : P ^ S and a = lo',a*A}g G T ^7r^^(^^* rxAdp)) 
the reduced section. Then 

d-^a + -a Aa = p*B, 

where p is the section of ttx^t. : S — X defined by p := 7rs,pO(T and B is the the Ad P-valued 
two-form induced on S by the curvature B = dA — {A A A) (z ^^{P, q)- 
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Proof. For U,V e (3.7) gives 

dMU. V) = V^(a(V)) - V^(a(f/)) - ai[U, V]), 

= la,d{a*A{V))U-[a*A{U),a*A{V)]l 

- ia,d{a*AiU))V - [a*AiV),a*A{U)]]^ 

-la,a*Ai[U,V])l 
= [a, d{a*Am V) - 2 [a*A{U), a*A{V)]h , 



a, a*dA - -a* {A A A) (U, V) 



la,a*B{U,V)l--{aAamV), 
p*B{U,V)-^-{aAamV). 



a, -a* (AAA) {U,V) 



Thus we obtain the formula 



d^a+^aAa = p*B 



as required. ■ 

Note that a G Q^{X,p* AdP), where p := tt^^p o cr and p* AdP is the pull-back vector 
bundle over X. Therefore rf-^a", a A a E Q"^ {X , p* Ad P) . Since the reduced curvature B 
belongs to the space fl^{T,,AdP), the pullback p*B e Q"^ {X , p* Ad P) . Thus, the formula 
d-^a + a A a = p*B is well-defined as an equality in f2^(X, p* Ad P). 

Corollary 3.2 Let he a smooth variation of the smooth section a : U ^ P and A a 
connection on the principal bundle vr^.p : P — S. Then 



De 



e=0 



crs{x),a;Al^ = V-^ la,A{6a)l^ + la{x),a*{is^B){x) 
-Hx),[A{6a{x)),a*A{x)]i, 



where B = dA — 1/2 A ^) is the curvature of the connection. 
Thus, the infinitesimal variations of a are of the form 



where p := tt^; L(^n* TXAdp) ° ^' '^^ '^'^ arbitrary variation of p vanishing on dU , r] is an 

arbitrary section in r(7rx,Adp) that projects to p and vanishes on dU , and B denotes the 
AdP -valued two-form induced on E by the curvature B. 

Proof. The second formula is a direct consequence of the first. To prove the first, one could 
verify the identity in local bundle charts. We prefer a global proof based on the previous 
lemma. 

Extending the bundle geometry in order to explicitly take account of variations achieves 
the objective. Consider P = M x P and X = M x X with the projection tt^ pi^iV) = 
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(e, 7rx,p(p))- Smooth sections of tt-^ p, a : X — P are in bijective correspondence with 
smooth variations of smooth sections of vr^ p, as follows: 

a{e,x) = {e,aeix)). 

Let G act on P by extending the action trivially to the M-factor. Thus, S:P/G' = ]RxS 
and AdP = M x Ad P. Since p := it j.o a it is clear that p{e,x) = {e,Ps{x)). Similarly, the 

connection A extends to a connection A G Q^{P,g) by setting A{de,Up) := A{up), for any 
Up e TpP. 

The section a : X — > P induces the reduced section a (see (3.2) for the general definition) 
whose explicit expression may be computed as follows: For {e^x) G X, m^, G T^-X, letting 
:= £a^, and using 

Tcr{de,u^) = {de,Tae{u^) + 5ae{x)), 

generates 

cr{de, Ucc) : = a{e, x),A iTa{de, u^)) = {{e, ae{x)),A {Tae{u^) + 5aeix))j 

IL JJ 

= (e,cx,{u,) + la,{x),A{6a,{x))}^^ . (3.11) 
The required formula is attained by evaluating the identity in Proposition 3.1, 

d-^d- + Ad = p*B, 

on the pair of vectors {0^,0), (dsyV^) for v^: G T^X. A direct computation shows that 
p*^m,0),{d„v^)) = (e, la,{x),l3{6a,{x),Tas{v,)l 



I {a A a) m, 0), {d„ V,)) = (e, la,{x), [A{6cx,{x)), A {Ta,{v,))]l 

To calculate d-^a {{8^,0), {ds,Vx)), let V G X(X) be such that V{x) = and use (3.9), 
(3.11), and [((a„0),(a„l^))] = Otoget 

d^am,0),{d,,v,)) = (v4,o) {Hde,V))) {e,x) - (V^,,^) (^(^,,0))) {e,x) 

-am,0),id,,V))] {e,x) 
= ^td,fi)i^,ae{V)){e,x) - Vfo,,,) [e, |(t,(x), ^((5(t,(x 

= -^{e,a,{v^)) - (^e, V;J [a,{x),A{5a,{x) 

The last three identities prove the first stated formula. ■ 

A covariant derivative on the tangent bundle te : TS — > S is needed in order to 
compute the variation of T^p. Given V^, (3.5) defines the V^-derivative which acts on 
functions X TT, and thus (3.6) provides the operator 

d^ := d^" : fii,,jx,rs) nl^Jx,Tj:) 
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defined by 

d^m V) = VBiKV)) - Vp(A(f/)) - A([f/, V]), (3.12) 
where A G fii^^(X,TS) and U,V e X{X). In particular, since J^S C L(7r^ ^^X, TS), 
sections of J^S ^ X are necessarily sections of Lln'^j^TXjT'E) X, tliat is, elements of 
f2^^^(X, TS) and thus operates on sections of the bundle J^S — > X. 
This differential operator satisfies the following property. 

Proposition 3.3 Let p : X ^ S be a smooth section of ttx,t.- Then 

where T'^{U, V) = V^V — Vyll — [U, V] is the torsion tensor of the connection V^. 

Proof. Recall that the section j^p : X — ^ J^S is interpreted in this formula in the following 
way. Given s G S, let x := Hx,T:{s) G X and so j^p{x) = T^p : T^X — >• Tp(a.)E, that is, one 
thinks of j^p as an element of f2^^^(X, TE). Given Ux,Vx G T^X and having chosen two 
vector fields U,V & X(X) such that U{x) = Ux and V{x) = Ux, (3.12) and (3.5) confer 



(f{j^p){Ux,Vx 



^BiTpiV)) 



V^iTpiU))-Tp{[U,V]) 



Dt 



Tp{v{c,m 



Dt 



Tp{U{c^m-Tp{[U,V]) 



where Ci(t) and C2(t) are curves in in X such that Ci(0) = C2(0) = x and Ci(0) = Ux, C2(0) = 
Wi;. Since p is a section of ttx.s, it is an embedding and the image p(X) is a submanifold 
of E. Thus, there exists vector fields U,V E X(E) such that U{p{x)) = Tp(U{x)) and 
V{p{x)) = Tp(y{x)). Accordingly, 

p*T''{ux,Vx) =T''{Tp{U{x)),Tp{V{x))) =T''{U,V){p{x)) 



[V^V-VfU-lU,?]) {p{x)) 



Dt 



mcim- 



D' 



t=o 



Dt 



Tp{V{c,{t))) 



Dt 
D"" 



u{p{c,m-[u.v]{p{^)) 



i=0 



Dt 



t=0 



Tp{u{c2m-Tpmv]) 



t=o 



which proves the statement. ■ 

The next result may be obtained using the previous formula by extending the bundle 
geometry to P as was done in the proof of Corollary 3.2 using Proposition 3.1. This time, 
however, we provide a different proof based on a standard formula for the torsion. 



Corollary 3.4 Let p^ be a smooth variation of the section p : U 
be a covariant derivative on TE. Then 



T,, Ux E TxX , and let 



D' 



De 



TxPeM = ^uJP + {Sp{x),Txp{Ux)) 



e=0 



where is the -derivative and T^{U,V) 
of the connection . 



[U, V] is the torsion tensor 
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Proof. This is a direct consequence of the formula 



where {e,t) E I x J a{e,t) G S is a smooth smooth function. Here it suffices to choose 
a{e,t) = Pe{c(t)), where c is a smooth curve in X such that c(0) = u^- ■ 

For simphcity, we will always choose a torsion free connection V^. In this case, the 
previous formulas simplify to 



d^{j^p) = and 



De 



£ = 



3.2 The Lagrange-Poincare field equations 

Let I: J^J:®^L {tx\.^TX, kdP) ^ M be the reduced La grangian (see (3.3)). This Section 
computes the Lagrange-Poincare equations given by the variational principle 







d 
de 



d 
de 



--qJu 



An affine connection on the affine bundle J^S ©s L {t:*x ^TX, Ad P) — S is required in order 
to obtain explicit formulas. Since the principal connection A brings a covariant derivative 
on AdP S, it suffices to choose a covariant derivative on the vector bundle TS S. 
This induces a connection on J^S given by 



V^^7:=ver^oV|7, 



(3.13) 



where ver''' is the vertical projection associated to the section 7 G r(7rs,jis), interpreted as a 
connection on tt^.e, and Z is a vector field on S. Here V'^ denotes the covariant derivative 
induced on L (^ttx y;TX,TT.) , from and from a covariant derivative V"^ on TX, that is, 



(vf£) (f/)(s) := vfmms) - i {VT.,Mas))U) 



(3.14) 



where £ is a section of the vector bundle L {n^^.TX^TT) ^ S, ^ G X(S), U G X(X), s G S, 
and recall that l{U){s) := ^(s)(t/(7rx,s(s))). However, the final result only depends on 
and not on V"^, see Janyska & Modugno [1996]. In this paper it is also shown that if is 
projectable onto a covariant derivative on X, then ^7 is an affine connection. 

Thus assuming a projectable covariant derivative is given, an affine connection on 
J^S L (tt^^^ 5^TX, Ad P) is obtained. Given a smooth function / on this affine bundle, 
define the fiber derivatives 



— ( 
5j^p 



d 

de 



lifpix) +ev,a{x)), 



e=0 



61 



[XI,W 



d 

de 



l{j^p{x), a{x) + ew) 



£=0 



where v G LiTr^X, \/p(2,)S) and w G LiT^X, AdPp(^x)) are arbitrary vectors. Note that 6l/6j^p 
and SI /5a are sections of the bundles L {■k\ j.T*X, VT.*) ^ X and L [ii\ j,T*X, AdP*) 
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X; both project to p. The derivative with respect to p is the horizontal derivative defined 
at {j^p{x),a{x)) by 

Jp^''^''') Te ' (3-15) 

where a{e) is a curve in S such that (j(0) = n G Tp(a.)S, and 5.(3,)) is the unique 

horizontal curve starting at (j^p(x), d"(a;)) and projecting to (J^e). 

Consider a variation ag of a given local section a : U ^ P and the reduced section (7^. 
Employing the decomposition of the e-derivative into its vertical and horizontal parts yields 



d 

de 



=oJu 



61 



K3^Pe{.x),ae{x))p 

, . d 



+ 



bp^ ''' de 
bp 



e=0 

' De 



dl 



U 



d 

de 

51 
6j^p 



{j^Peix),a,{x)) ] p 



e=0 



De 



J^Peix) ) p 



(3.16) 



e=0 



ae^x) ) p 



£ = 



where D-^^^/De and D^ /De denote the covariant derivatives associated to the connection 
V"' ^ on J^S S and to the induced covariant derivative on L {ji^jJ^X^KdP^ S, 
respectively. 

The second term may be computed using the following relation: 



D-^ 



De 



D 



J Pe[X) [V^j 



De 



e=0 



{fpe{x){v^)) . 



(3.17) 



e=0 



This relation is obtained from the definition of the induced covariant derivative D^ /De on 
L{t[\jJ^X,TY?). Given a curve 7^ G J^S, (3.14) shows that 



D^ 



e=0 



Is {Vx) 



D 



TS 



De 



e=0 



{le-Ve) -70- 



D 



TX 



De 



e=0 



where Vg G TX is a curve such that G T^^X and G X is such that 7^ G LiT^^X, T^^S). 
In the present case 7^ = j^ps{x) and variations in TX are not considered, so = x and 
Ve = Vx- Thus 

U'Psix)) (Vx) 



[iDe 



e=0 



De 



{fps{x)-Vx). 



(3.18) 



£=0 



Denoting the connector map of by K^^ and recalling that is projectable allows the 
following calculation: 



Tnx,s 



D 



TT. 



De 



ij^Peix)-V^ 



e=0 

TTtix,t. 

d 

T ^ 



Ttt 



d 

de 



ij^Peix)-v^ 



e=0 

0. 



d 

de 



ifpeix)-Va: 



e=0 
d 

de 



Tirx,T {j^Peix)-v^) 



e=0 
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This proves that the expression (3.18) is vertical. Thus, by the definition (3.13) of V"^^^, the 
identity (3.17) is proved. 

The third term in equation (3.16) may be evaluated using the equality 



[iDe 



(TAX [V, 



e=0 



(3.19) 



£ = 



whose proof is similar to that of (3.17). 

Using (3.17), (3.19), and Lemmas 3.2, 3.4, the expression (3.16) may be rewritten 



d 

de 



=oJu 
51 



51 



5p' 



u 



+ 



51 



^-,V^r/-[r/,a]+p*(i,,S))/x 



51 _ ^.^s _/ 
5p 5j^p \5a 



51 



51 



+ 



Since 5p and t] are arbitrary, this results in the vertical and horizontal Lagrange-Poincare 
equations given by 



. A SI 51 . SI ^ ^ 51 / 51 . 

div-^ — - ad; — = and div^ — + ( — , i^pT^ 

OCT 5a 5p 5j^p \oj p 



51 
5^ 



/iTpB), (3.20) 



respectively, where the second equation has to be considered as an equation in V^(2,)S*. For 
simplicity, we will suppose that is torsion free. 

Here div^ denotes the divergence associated with V"^, 



div-^ : r ( TT, 



which is defined as minus the adjoint differential operator to V"^: 



X 



div^/(x),^(a;))/i 



X 



for all / e r (vr^,i(,^^T*x,Adp.)) and ^ e r(7rx,Adp) such that vr^,i(,^^r*x,Adp.) ° ^ 
TTs.AdP o C,- In the vertical equation, ad* denotes the map 

-> r (vTx.AdP*) ) At) ^ ad; /i. 



ad* : r f vr 



'X,L(7rJ,_5,TX,AdP 

well-defined when 



X r vr 



j:,L{Tr* ^T*X,Ad P* 



O p = TT. 



S,L(7r* ^TX,AdP 



o a. Similarly, the operator 



div^ : r TT 



(3.21) 
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is the divergence associated to the V^-derivative restricted to vertical valued sections: 

: r(7rx,ys) ^ T (7r^,i(^i ^.Txys)) ■ 

Note that such a restriction is possible since is projectable. The results obtained above 
are summarized in the following theorem. 

Theorem 3.5 Let 7Tx,p : P ^ X be a locally trivial fiber bundle over an oriented manifold 
with volume form fi. Let L : J^P —>■ be a Lagrangian which is invariant under a free and 
proper left action ^ : G x P —>■ P such that 

T^x,p °^g = T^x,p, for all g EG. 

Let 7rs,p : P — S := P/G be the associated principal bundle. 

Fix a connection A on t^y.,p (ind let I : J^T, ©s L (tt^ ^^TX, Ad P) ^ M. be the reduced 
Lagrangian induced on the quotient by means of the identification (2.8). Let a : U P be 
a smooth local section of tix,p, define the reduced local section a ofn-^^^^, TXAdp) 

aix) = laix),AiTM-)],, 

and the local section p := 7rs,p o a of -nx,Ji- Fix a projectable covariant derivative on TT, 
and suppose, for simplicity, that is torsion free. Then the following are equivalent: 

i The variational principle 

S [ L(j V)^ = 0, 
Ju 

holds for arbitrary vertical variations 6a vanishing on dU . 
ii The section a satisfies the covariant Euler- Lagrange equations for Lfi. 
ii The variational principle 

Jx 

holds, for variations of the form 6-^a = V-^f]—[f], a]+B{6p, Tp), where 6p is an arbitrary 
variation of p vanishing on dU and fj is an arbitrary section of it x, Ad p vanishing on 
dU and such that vr^^AdP o fj = p. 

iv The section a satisfies the Lagrange-Poincare field equations 

61 ^ 61 161 . ~\ 
— -div -— = ( j^,iTpB ) , 

op 6fp \6a I ^2 22) 

div-^ 1^ - ad; 1^ = 0. 
6a 6a 

In the case of a connection with torsion, a term involving the torsion tensor has to 
be added in the horizontal Lagrange-Poincare field equations, see (3.20). 

The Lagrange-Poincare field equations are now examined in the particular cases men- 
tioned before. 
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If G = {1} then H = P, p = a, I = L, and there is no reduction. In this case (3.22) 
becomes 

— - div^ = 0, 3.23 

which is just a restatement of the covariant Euler-Lagrange equations, using a pro- 
jectable and torsion free covariant derivative on TP P. 

If 7Tx,p : P ^ X is a principal bundle and the symmetry group is the structure group, 
then E = X and the section p is absent since it is the identity on X. Therefore, the 
reduced variation reads = V'^fj — [fj, a], where fj is an arbitrary section of 7rx,AdP5 
and the Lagrange-Poincare field equations (3.22) read 

div-^ — - ad; — = 0. 
da da 

Thus the covariant Euler-Poincare equations are recovered; see Theorem 3.1 of Cas- 
trillon-Lopcz ct al. [2000]. 

If 7rx,p : P ^ X is a principal bundle whose structure group contains the symmetry 
group as a subgroup, the equations (3.22) coincide with the equations (4.11) obtained 
in Castrillon-Lopez & Ratiu [2( 



If P = M X Q where t^m^q : Q M is a G-principal bundle then S = M x M. The 
sections a G r(7r]R,p) and p G r(7rK,s) read a{t) = {t,q{t)) and p{t) = {t,m{t)), where 
m(t) = 'n'M,Q{(l{t)) G M. The first jet extensions j V and j^p are identified with {t, q{t)) 
and (t, m(t)). 

In this particular situation, the connection ^ on P is always chosen to be induced by a 
connection 7 on Q. In this case, the reduced section a is identified with [t, v{t)), where 
v{t) = [q{t),-f{q{t))]c. Similarly, a section ^ G r(7rM,Adp) covering p{t) = (t,m{t)) 
reads ^{t) = {t,^{t)), where ^(t) G Ad^(t)Q- The V^-derivative of ^ can thus be 
identified with the covariant time derivative (D'^/Dt) ^{t). Using all these observations, 
the second equation of (3.22) reads 

and the variation of v is = — [f],v] + B{6m,m), where B is the reduced 
curvature of 7. Recall that writing the horizontal equation requires a projectable 
covariant derivative on TS, which is also assumed to be torsion free for simplicity. 
In this classical case, the covariant derivative is constructed from a torsion free covariant 
derivative V on TM and the natural covariant derivative on TM. In this case, is 
obviously projectable and torsion free. The first equation of (3.22) reads 

61 61 _ / 61 ^ ~ 

6m Dt 6m \6v^ ^ 

Thus the classical Lagrange-Poincare equations obtained by standard Lagrangian re- 
duction are recovered; see Theorem 3.4.1 in Cendra et al. [2001]. Note that here the 
Lagrangian is allowed to be time-dependent. 
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If V has a torsion T^, the horizontal equation reads 



61 61 / 61 . v\ /^^ ■ « 



6m Dt 6m \ 6rn I \6v 

see (3.20). 

In the particular case, G = {1}, there is no reduction and the vertical equation is 
absent. In this case the horizontal equation reads 

6L 6L 



Dt 6q 6q 



0. 



Of course, this recovers the standard Euler-Lagrange equation written with the help 
of a connection. In the case when the connection has torsion, this reads 

6L 6L / 61 . _v 



Dt 6q 6q \6q 

see (3) in Gamboa & Solomin [2003]. Recall that the usual way to write the Euler- 
Lagrange equations 

ddL dL _^ 
dt dq dq 

makes sense only locally; see (3.1). 

Another particular case arises when Q = G. In this case, there is no horizontal equation 
and the vertical equation gives the Euler-Poincare equation. Indeed, in this case, all 
the connections are equivalent (the bundle Q ^ Q/G is over a point) and the covariant 
time derivative on the adjoint bundle becomes the ordinary time derivative on the Lie 
algebra q. These observations and (3.22) recover the Euler-Poincare equation together 
with the associated constrained variations 



3.3 Reconstruction 

Having derived the Lagrange-Poincare field equations it is natural to turn to the problem 
of reconstruction of solutions to the original Euler-Lagrange equation from solutions to the 
reduced equations. More precisely, given a solution section a : tJ ^ L{T:\j.TX,KdP) of 
the Lagrange-Poincare equations, how can one construct a solution section a : If P 
of the Euler-Lagrange equations? Note that Theorem 3.5 does not consider this problem, 
since the section a is given a priori. This section deals with the reconstruction problem and 
demonstrates that reconstruction of field theories requires an extra integrability condition. 



Induced connection. A section a G ^{'n'x,L{TX,AdP)) induces a G- principal bundle 
X and a connection A"^ on it as follows: The subset P^ G P is defined by 

P^ := n^UpiX)) = {peP I 7rs,p(p) G p{X)} , 
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where p := tTj, ^j^^* tx Adp) ° ^ r(7rx,s)- Since p is a section, it is an injective immersion 
and a homeomorphism onto its image. Thus, the image p{X) is a submanifold of S. Now, 
since tt^^p is a submersion, it is transversal to the submanifold p{X). This proves that is 
a submanifold of P, whose tangent space at p is 

TpP^ = {Tp7rj,,p)-' (r,p(T,S)), s = 7rs,p(p). 

The manifold may be endowed with the structure of a G-principal bundle over X by 
restriction of the G-action on P. Note that can be identified with the pull-back bundle 
p*P = {{x,p) I vrs,p(p) = p{x)}, the identification being given by 

peP"^ (7rx,p(p),p) G p*P 

The section a may be regarded as a section of the vector bundle L{TX, p* Ad P) ~ L{TX, Ad P''), 
and thus induces an equivariant and vertical one-form u"' G Q^{P^,q). The isomorphism 
ni(P'',0) ^ L(TX, AdP^) is written explicitly as follows: 

LO^ e iF(P^g) ^ a G L(TX, AdPO, := b,^;'^^)!^, (3.24) 

where Ux = T7Tx,p{up) G T^X, Up G TpP. The connection A on tt^^p : P ^ S naturally 
induces a connection on P^ — > X. A new connection A" is thereby obtained on P^ ^ X. 
Concretely, 

A^ := A" - cu^. 

Thus one may interpret the vertical solution of the Lagrange-Poincare field equations as 
describing an affine modification to the a priori connection A''. The modified connection A'^ 
is the correct choice of connection for reconstruction, as is explained below. 

Reconstruction condition. We now prove that if a is the reduced section associated to 
a section a G r(7rx,p) then A'' is flat. Indeed, in this case P^ = {$g(o"(x)) \ g E G,x E X} 
and for p = cr{x) G P^ and Vp G TpP^ formula (3.24) gives 

A^'ivp) = AP{vp) - u^{vp) = A{vp) - AiT,a{Tp7ixAvp))) (3-25) 

since a{x) = |cr(x), (T*^(x)]g for all x E X. Recall that Up G TpP^ if and only if Tp'UY.^piUf^ G 
TxpiTxX). That is, in terms of cr, 

TpHj^^piup) G T^{x)'^s,p {TxaiTxX)) . 

This proves that T^(x)P^ = TxCr(TxX) + V^(^x)P at p = cr(x), where V„(^x)P = ker(T^(^)7rs,p) 
is the vertical space relative to 7rs,p. Thus, for p = a{x), any Vp G TpP^ reads Vp = 
Tx(j(vx) + ^p(p)- Inserting this expression for Vp into (3.25), reveals the condition ^ = 0. 
This proves that the ^"^-horizontal subspace at a{x) is given by 

Hf;)P = TxaiTxX). 

This horizontal distribution is integrable, the integral leaves being given by {^g{a{x)) \ x G 
X} = $g(Im(cr)), for each g E G. Thus, the connection A" on P^ is fiat and the horizontality 
condition 

a*A'' = 0. (3.26) 
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is a necessary condition for reconstruction. 

Conversely, consider a section a of x,l{tx, Ad p) such that the connection on P'' is 
flat and has trivial holonomy. Since the connection A'^ is flat, the horizontal distribution is 
integrable and the leaves cover the base, that is, given a leaf L, each fiber intersects the leaf 
L at least once. Since the holonomy is trivial, each fiber intersects the leaf exactly once. This 
construction shows that each integral leaf of the horizontal distribution defines a section of 
the bundle X. Thus a family of sections $c, o a of itx,p that project via n^.^p to p is 

attained. Since 

la, a*Aj^ = la, a* A' + uj% = {a, u% = a, 

the section a is the reduced section associated to the family of sections ^goa for each g E G. 
The horizontahty condition (3.26) is, of course, satisfied. 

Recall that the flatness of the connection does not imply that the holonomy is trivial 
unless the base is simply connected or the holonomy group is connected. Note that this fact 
implies that the holonomy of a flat connection is locally trivial, that is, for every x E X, 
there exists an open neighborhood U such that the holonomy of P\u is trivial. 

The situation is summarized in the following reconstruction theorem. 

Theorem 3.6 Fix a connection A on the principal bundle n^^p : P S, consider a G- 
invariant Lagrangian L and the reduced Lagrangian I. 

If a : If ^ P is a solution of the Euler- Lagrange field equations, then the reduced section 
a is a solution of the Lagrange- Poincare field equations. Moreover the connection A'^ on P'' 
is flat and the horizontality condition (3.26) holds. 

Conversely, given a solution a of the Lagrange- Poincare equations on U such that A" is 
flat and has trivial holonomy over an open set containing U, the family o a , g E G, of 
solutions of the Euler-Lagrange field equations are given by the integral leaves of the horizontal 
distribution associated to A" . In addition, the horizontality condition (3.26) holds. If the 
connection A" is flat one can always restrict it to an open simply connected set contained in 
U so that its holonomy on U is automatically zero. 

Note that the curvature of A'^ is B — d-^u" — (1/2) u;'^ Acj"". Therefore, the reconstruction 
condition is 

B - d-^a;'^ -^u"" Au"" = on P". (3.27) 

This condition has to be seen as an equality in the space Q'^{P'',q) of equivariant vertical 
two- forms. The isomorphism (3.24) shows it is equivalent to assume that the corresponding 
two-form in n'^{X,AdPP) = n'^{X,p* AdP) vanishes. Applying (3.24) to equation (3.27) 
recovers the formula 

d-^a +^a Aa = p*B. (3.28) 

Reconstruction equation. When reconstructing solutions of the Euler-Lagrange field 
equations one needs to add (3.28) to the reduced field equations (3.22) since there could 
be solutions to the Lagrange-Poincare field equations (3.22) that do not correspond to the 
original Euler-Lagrange system. Given a solution (p, a) as specified above, (3.24) uniquely 
determines u'^ by the formula 

a{u,) = [p, {RoT^{T,p{u,)))] , p G n^^p (p(x)) , G T,X, (3.29) 
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since TpTix,p {}^oip{Txp{ux))) = Ux- Thus, uj" is completely determined in terms of {p,(t)- 
For a section a G F (ttx.pp), the horizontality condition (3.26) for is 

= = a*A^ - a*uj^ = a*AP - p* (Hor^)* cj'^ (3.30) 

because ixi" (T^.a(M^.)) = (Hor^^^) (T^p(m^))) since G fii(P^,0). Note that following the 
determination of io" by (3.29) the only unknown quantity in (3.30) is a. We now show that 
(3.30) gives a first order PDF that determines a. 

If Ux G T^X, then by (3.30) and the horizontal-vertical decomposition relative to the 
connection A'^ ^ 

Txcr{ux) = Hor^(''^) (T^(x)7rs,p (TxCr{ux))) + {A^ iTxCr{ux)))p{cr{x)) 
= Hor^f,) (Txpiux)) + ia*APiux))p (a(a;)) 
= Hor^f,) (Txpiux)) + (p* (Horf )* ^-)^ (a(x)). 

This gives the following first order reconstruction PDF for a: 

Txa = Hor^(,.) o T.p + (^^ o Hor^f,) o T.p)^ (a(x)). (3.31) 

Theorem 3.6 can now be interpreted as asserting that given a solution (p, a) of equations 
(3.22) and (3.28), there exists a unique solution a to the reconstruction equation (3.31) in 
a neighborhood where A'^ has trivial holonomy. This section a solves the corresponding 
Fuler-Lagrange equations for the unreduced problem. 

As a final comment, note that (3.31) is the field theoretic analogue of the classical recon- 
struction equation ^ = g^^g associated to the Fuler-Poincare equations. 

Particular cases. The reconstruction condition specializes to the particular CclSGS clS fol- 
lows: 

i It G = {1} there is no reduction and, therefore, no reconstruction condition. 

ii In this case the variable p is absent, so P'' = P. Moreover, the reduced section a 
turns out to be associated, via the map to a section C of J^P/G X, that can be 
interpreted as a connection on P. This connection ( does not depend on the chosen A 
and A'^ turns out to be the connection one-form associated to (. The reconstruction 
condition is simply that the curvature of this connection (or of () is zero. This recovers 
the reconstruction condition that in the case of covariant Fuler-Poincare reduction; see 
§3.2 of Castrillon-Lopez et al. [2000]. 

iii The reconstruction condition is the same as in Castrillon-Lopez & Ratiu [2003]. 

iv Since X = M, the base is one-dimensional and every connection is flat. Since M 
is simply connected the holonomy is trivial. The reconstruction condition is always 
satisfied. This agree with the fact that in classical Lagrangian reduction, the solution 
of the Fuler-Lagrange equations can always be constructed from that of the reduced 
equations. 
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4 Conservation laws and representations 



In applications there is often a natural choice of gauge that is used to formulate the Lagrange- 
Poincare field equations in a convenient local form. 

This section describes the two predominant choices of representations that occur, the 
spatial representation and the convective representation. The Lagrange-Poincare equations 
(3.22) are given locally using these choices of gauge. The spatial representation yields 
Noether's Theorem as the vertical equation whilst the convective representation has the 
Euler-Poincare equation as its vertical equation. This observation shows that the Lagrange- 
Poincare equations are equivalent to Noether's Theorem, a statement often found in the 
literature when dealing with concrete applications. 

This section also formulates a global version of the Kelvin-Noether Theorem that gener- 
alizes the result for classical systems given, for example, in Ccndra et al. [1998]; Holm et al. 
[1998]. 



4.1 Representations and Noether's Theorem 

A section a eT (tt^.p) introduces a representation of Ad* P which, in turn, yields local equa- 
tions for the vertical part of (3.22). The two natural choices of section and their associated 
representations are described below. 



Convective representation. Suppose one seeks a local solution a : U C X — P of (3.22) 
and (3.31) in a trivialization of P over U <Z X. Let p := vr^^p o a and V := p{U) C S. 
Suppose further that a flat connection A exists o\i P ^ V . Then, there exists a unique 
section 7 : V P such that TsjIvs) € HsP, for all s G S and Vg G T^S. Therefore, the 
section o"^ := 7 o p g T^nup) has the property that T^cr^ (v^) G H^hf^^^P for all G TX\u. 
Such a section is called a horizontal section. 



Remark 4.1 It may not be possible to find a flat connection A on an arbitrary open set 
V <Z Tj. The convective representation is not defined in such cases. In applications one 
may find that shrinking the set C X yields a suitable V <Z Tj such that the convective 
representation makes sense. 

The map g : U C X ^ G such that a{x) = ^g(^x)0'^{x) for all x E U together with (2.4) 
produce 

Xa^ (^) = Xa^ [k, ^*^lg) = g'^dg =:^en^ (X, 0) . 
Consequently, the vertical part of equations (3.22) composed with x^-h yields 



Thus, the local representation of the vertical Lagrange-Poincare equation in this gauge is 

div|-ad-|^0 (4.1) 

which recovers the Euler-Poincare equation. This choice of gauge is called the convective 
representation, see Ccndra et al. [1998]. 
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Spatial representation. With the same notation as for the convective representation, 
a = ga^, the map Xa applied to a yields 

X. (^) = Xa = dgg-^ =: E e (X, q) . 

Accordingly, 




Thus, the local representation of the vertical Lagrange-Poincare equation (3.22) in this gauge 
reads 

div-^ = (4.2) 

0^ 

which is Noether's Theorem. This choice of gauge is called the spatial representation. Note 
that Adg^ = S. Therefore both the Euler-Poincare equations and Noether's Theorem are 
local representations of equations (3.22) corresponding to a particular choice of gauge. In 
particular, the Euler-Poincare equation is equivalent to Noether's Theorem. 

Remark 4.2 When the convective representation can not be defined it is still possible to fix 
S = a* A and proceed with the construction of the spatial representation without the use of 
a^. Thus the spatial representation is always well-defined, while the convective representation 
is not. 

Remark 4.3 In classical Lagrangian reduction when X = M it is always possible to con- 
struct a local horizontal section. Therefore the convective representation is always well- 
defined for classical systems. 

4.2 The Kelvin-Noether theorem 

Given any manifold C on which G acts, the associated bundle is a fiber bundle over S defined 
by 

He ■.= PxGC = iPxC)/G, 

where the action of G on P x C is the diagonal action. The adjoint and coadjoint bundles. 
Ad P and Ad* P are associated bundles with C = q and C = q* respectively. The action of 
G on for AdP is the adjoint action whilst the action on G on q* is the coadjoint action. 
The equivalence class of (p, c) G P x C is denoted 

|p,c]cGPxgC. 

The lifted action of Gon TC enables the definition of Tire = P TC. The infinitesimal 
action Ad P x He — > Ti-rc on Tic as follows: 

b,ei,-b,clc = b,ec(*c' (4-3) 

where the vector field ^ ^(C) denotes the infinitesimal generator of ^ G g on C. 
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A connection form A on tt-e^p yields a covariant tangent functor T-^ defined on sections 
of TixMc by 



TC ' 



(4.4) 



Therefore, if c is a section of iTx,nc covering the section p of nx,j:, then T-^c{vx) G {Ji-Tc)p{x)- 
The map K can be described by a G-equivariant map K. : C ^ q** defined by the relation 

ir(b,cy = b,/c(c)i^„. 

Here Q** denotes the double dual of the Lie algebra. For an example where the distinction 
between g and g** arises, see Holm et al. [1998]. The derivative of K : Tic Ad** Q may 
be defined as follows: 

Note that dK : (7^tc)s (Ad** P), and d^ {K o c) = dK o T^c, where c = {a, cj^ denotes 
a section of 'n'x,Hc- Indeed, 

d-^ {Koc) = [a, d (/C o c) - (a*^) ■ (/C o c)]g,, 
= |cr, die oTc — die o (cr*^)|^ o c]^,^ 
= [a,rf/Co(rc-(a*^)cOc)]^.. 
= rf/sT o T-^c. 

This relation is described by the following commutative diagram: 

r (tIx,L{-^x,sTX,Htc)) r (nx,L{n*^^TX,Ad" P)) 



^A 



TiTlX 



He 



K 



TiTlX 



Ad" P 



The infinitesimal actions described in (4.3) on Ad** P and Tic are readily observed to be 
related via 

lp,a,-K{c) = dK(lp,n.c). (4.5) 



(4.6) 



Additionally, observe the following relationship: 

div (z/ , /i) = (^d'^u , /i) + (z/ , div^/i). 



where u and p, are sections of vrx,Ad** p and tixm-^*^ ^T'^AdP*) respectively, and both cover 
the same section p of 7rx,s- 

The Kelvin-Noether Theorem may be stated as follows: 

Theorem 4.4 Let (p, o") G F (vrx,s) x (X, AdP) be a solution to the Lagrange-Poincare 
equations (3.22), and c G F (ttx^Hc) cover p while satisfying 



T^c + ac o c = 0. 



(4.7) 
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// K : Tic Ad**P fiber-preserving map that covers the identity on S then the associated 
circulation 

I := (koc, eX(X) 

satisfies 

divJ = 0. (4.8) 

Proof. The result is obtained via a direct calculation that uses (3.22), (4.6), (4.5), and (4.7) 
as follows: 

div / = div ( K o c 



.-(Koe),|) + (Koe,d.v-|) 
{d^ + adr) {K o c) ,^) + (ko-c, (div^ - ad;) 



6ct/ \ ' ^ 5a/ 
{d^ + adr) (Koc),^ 

dK o (T-^c + ac o c) , — \ = 0, 
as required. ■ 

Recall that classical Lagrangian reduction (particular case iv), used the formulation 
ait) = (t,q(t)) and a{t) = lq{t), A{q{t))}^. In this case (4.7) becomes 

T^c + acoc=lq,c-{A (g))^ o cj^^ + |g, {A (g))^ o cj^^ = |g, c] = 0. 

Therefore (4.7) diminishes to the assumptions for the classical Kelvin-Noether theorem; see 
Holm et al. [1998]. Furthermore the conclusion to Theorem 4.4 in this context becomes 

These results extend those of Cendra et al. [1998] to the Lagrange-Poincare context. 



5 Applications 

This Section presents a brief outline of some applications of the Lagrange-Poincare field 
equations. 

The first application is the minimal immersion problem, which is treated explicitly in 
local coordinate form. 

The second and third applications are classical and covariant metamorphosis image dy- 
namics treated in coordinate-free form. The classical formulation may be understood as a 
change of variables from the treatment given in Holm et al. [2008]. This coordinate transfor- 
mation formally decouples the equations. The covariant formulation provides a new insight 
into the problem. This covariant formulation aims to provide a basis for the future use of 
multisymplectic integrators in image dynamics. 
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5.1 Minimal immersions 



An interesting classical problem that has applications ranging from the shape of soap bubbles 
to string theory is that of minimal embeddings or, more generally, minimal immersions. 
In string theory the Nambu action describes a world sheet in spacetime; see for example 
Nakahara [2003] . The soap bubble problem is a generalization of the isoperimetric problem 
studied by Newton amongst others. 

Given a manifold X and a pseudo-Riemannian manifold {Q,g), the problem is to find an 
immersion rj G Imm(X, Q) such that the surface area of rj {X) is minimized. For the soap 
bubble problem there is an additional constraint on the volume enclosed by t] (X) which we 
shall not treat here. 

The minimal immersion problem may be cast into the bundle picture described in §2 as 
follows: Let P = X xQ and vrx,p = Pi be projection on the first factor. Since vrx,p is a trivial 
fiber bundle, sections a E T (vrx,p) of nx^p can be represented by smooth maps r] : X ^ Q, 
namely a{x) = {x, r]{x)). Further, consider sections a G F {ttx,p) such that rj G Imm {X, Q) 
is an immersion, that is, a G F^"*™' (7rx,p) where 

F^"'" (vrx,p) = {a G F (7rx,p) h G Imm (X, Q)} . 

Since rj is an immersion, h = ri*g is a metric on X. Locally h reads 

hij = {gafs o r])dir]°'djr]'^, 

where the indices ... denotes coordinates on X and a,j3, ... are coordinates on Q. The 
Lagrangian density for the minimal immersions is the volume form on X associated with the 
metric h = ri*g: 

C:J^P^ A"+iX, C (j V) = ^{d'^+^x =: L{fa)d''+^x, 
where \h\ := | det h\. Locally, the Lagrangian reads 



L (xi,qa,i^j^ = det (^g^p{q)u°'u^- 



where (xj, q^, i^j) are the natural coordinates on J^P. 

Since the bundle P is trivial, the covariant Euler- Lagrange equations read 

^ - div*^ ^ = 0' (5-1) 
07] oj'-r] 

where div*^ is the divergence operator associated to a V*^-derivative (see (3.5)). Since Q 
is a pseudo-Riemannian manifold, the Levi-Civita connection provides a natural choice of 
covariant derivative on TQ. The covariant Euler-Lagrange equations may be calculated by 
use of the following formula for the derivative of the determinant of an invertible matrix 

D det{K) ■ 6K = (det K) tr {K~^5K) . (5.2) 

Since V*^ is the Levi-Civita covariant derivative, the first term of (5.1) vanishes and the 
covariant Euler-Lagrange equations have the local representation 

Vf = 0, (5.3) 
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where 

8L 1 
Pa = ^ = sign{h)\h\^h'^ g^pdjT]'^ 

and V*^ is the V'^-derivative. 

Note that when X = M then equation (5.3) reduces to 

which is just the geodesic equation on {Q,g). 

Now consider the case when the isometry group G = Iso(Q, g) of g acts freely and properly 
on Q. Then 71m,q : Q M := Q/G is a. principal G-bundle. Whereupon the group action 
G X P — >• P by (/, (x, q)) (x, fq) gives P a principal G-bundle structure over S = X x M. 
The geometric setup is as described in §2.1; this fact is elucidated by the following diagram: 

X xQ ^ X 




X X M 



The diagram also reveals that reduction of the minimal immersion problem is an example 
of fiber bundle reduction, the natural extension of the classical Lagrange-Poincare reduction 
discussed in §1.3. 

Identification p = vrs^poa G F^™™ (^x,e) with r G Imm (X, M) may be effected by writing 
p{x) = {x,r{x)), where r = tcm^q ° V- 

Lemma 5.1 The Lagrangian density C : J^P A""^^X defined by 

is G-invariant. Fixing a particular principal connection A on P, the reduced Lagrangian 
I : {J^P)/G —>■ M may be expressed as 

^ (jV>o-) = h = r*gM © ^*gAdQ, 

where r := 'Km,q ° p, Qm is the Riemannian metric on M defined by gM = (Hor"^)* (7, and 
Qaaq is the vector bundle metric on AdQ defined by 

QAdQim) (\q, , |g, := g{q) {^Q{q),VQ{q)) ■ 
Proof. Since g is G invariant in the sense that f*g = g for all f E G, 

^ (/ ■ jV) = I {fvYgl'^d^+'x = \r^*rg\'^d-+'x 

= \rfg\^d''-^^x = C{3^a). (5.4) 

Thus C: J^P ^ A"+iX is left G-invariant. 
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Since Q is a pseudo-Riemannian manifold there exists a mechanical connections whose 
horizontal spaces are orthogonal complements of the vertical spaces. This connection on Q 
induces a unique connection A on P. The reduced configuration space J^P/G is identified 
with J^S©L (vr^s^ 5^TX, Ad P) using the isomorphism /?_4 defined in (2.8). Consequently, the 
reduced Lagrangian / : J^P/G = J^T, ©s L (^n'^ j.TXjAdP) M may be expressed by 

I {Pa [7]) d^'^^ = L (7) , for all 7 G J^P. (5.5) 

Given a section a{x) = {x, rj{x)), the objective is to compute h := rfg in terms of the reduced 
quantities p(x) = (x,r(x)) and a{x) = la{x),a*A{x)}^ = lr]{x),ri*A{x)}^, where in the last 
equality, the adjoint bundles of P and Q have been identified in the canonical way. Since A 
is the mechanical connection, 



9{(l){uq,Vg) = gM{Tr{q)){TTT{uq),TqTT{vq)) + gAdQ{Tr{q)) {q, A{ug)j Aq, A{v. 



Thus, 



r]*g{x){u^,Vx) = g{r]{x))(T^r]{ux) ,T^r]{v^)) 
= gM{r{x)){T^r{u^),T^r{v^)) 

+ gAdQ{r{x)) ([t]{x),A {T,^r]xM)\^ , {T^vivx))\ 
= ir*9M) (mx, Vx) + 9AdQir{x)) {a{u^), a{v^)) , 
which is more compactly expressed by 

V*9 = r*gM © ^*9AdQ- 
This formula, together with (5.4) and (5.5) combine to show that 

/(jV>^) = l^l^ h = r*gM®^*gAd 
which is precisely the statement that was to be proved. ■ 

Setting := V"^ © V^^, where V^^ is the Levi-Civita connection associated to gM 
on M = Q/G and V"^ is an affine connection on X, formula (5.2) gives the functional 
derivatives, 



61 



d 
de 
1 



e=0 



sign{h) y' \ h\ tr I h 



d 
de 



guir) {Tr + ev, Tr + ev] 



e=0 



and 



51 



,w 



Slg 

d 

de 



n{h) ^/\h\ tr {{r*gM) ^ guir) {Tr,v)) 



I (jV,o- + ew) 
^sign(/i)A/|7z]"tr {h^^ 



de 

-1 



9AdQ{r) (0- + ew,a + ew] 



e=0 



sign(/i)vR^tr ((cr* 5(Ad) 9AdQ{r) (o-,^)) . 
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Derivation of the final functional derivative, 61/ 6p, requires a curve cr (e) G S. The horizontal 
curve 

then denotes the horizontal lift of a{e) with respect to the affine connection V"'^^ © V"^. The 
final functional derivative is then defined by the relation 



61 



d 



6p^V^''^ de 



£ = ^ ^ 



as in (3.15). Since V*^ is the Levi-Civita connection with respect to qm there is no contri- 
bution to 61 /6p from the qm terms of /. Thus, 



61 \ . ^ d 



1 



e=0 



d 



n 



e=0 



where 



7^ = tr 



cr, 



-1 r 



Oip(x),s(x))) 9AdQ{r){0) ((^iPpix),aix))) 9AdQir,){e) 



Since V^a^i^^^-) -^^-^-^ = by definition, 71 G cr*JF (S). Therefore introducing 71 E J-' (S) such 



Note that the first factor in the trace is evaluated at e = before the derivative is taken. 
Since V:^a^.i , ^ 

that a*7l = 71 gives 



d 
de 



(5.6) 



e=0 



This leads to the formula 
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u 



6p' ~ I 2 

All the functional derivatives collected together read 

61 
6(r 



sign(/i)V|/i|M Jc^ 7^. 



61 
6fp 



sign(/i)v^(a*(?AdQ) '[a,a*iIA)]a=:n 
sigI^{h)^/\h\{r*gM)~^ gu (r) (Tr, ■) =: P 



(5.7) 



61^ 
6p 



lsign{h)y^\d^7Z 
where I : Q x g ^ g* is defined by 

(I(g)^, T]) := giiqiq), r]Q{q)) =gAdp[ [q, , [g, rj\^ 



38 



Now, restricting attention to a local neighborhood U C X consider cr as a local section of 
the principal bundle over U. Using capital letters for g coordinates which are raised and 
lowered by I, lower case letters for X coordinates raised and lowered by h and greek letters 
for M coordinates raised and lowered by qm results in the following local representations: 



where 



61 
51 

51 
Tp 



C 



a, {a* A)' 



C {a* A)' =: 



Cr: =: PI 



2C 



C := sign(/i)\/|7i|. 



The Lagrange-Poincare equations read 

C 
" 2" 
adm 



div'^P 



div-^n 



0, 



(5i 



(5.9) 



where B is the reduced curvature form associated to the connection A. Equations (5.9) can 
be written in the local coordinates as 



^ I a 



2C 



(5.10) 



where B'^g is the local representation of the curvature form B associated to A. The right hand 



side of the first of equations (5.10) measures the deviation of r being a minimal immersion 
m. M = Q/G whilst the second of equations (5.10) is Noether's Theorem. 

For reconstruction, consider the form uj'^ defined by equation (3.29). The curvature 
relation (3.27) is given in coordinate-free form as 



B 



'-u^ A cj^ = 0, 



which is expressed in local coordinates by 
, ...J 1 



(to' 



la + '^KM 



,^ )a )p 



(5.11) 



where cj^^ are the structure constants for g. The Reconstruction Theorem 3.6 states that if 
one can solve equations (5.10) and (5.11) then there exists a unique solution of (3.31) which 
is also a solution of (5.3). 

The degree of geometric content in equations (5.9)-(5.11) indicates why the problem of 
minimal immersions has fascinated mathematicians and other scientists for so long. One may 
expect that further investigations of this problem will continue to produce rich mathematical 
results. 
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5.2 Metamorphosis image dynamics 

The metamorphosis framework is an interesting approach to the control theory problem of 
how best to match one image to another, particularly when the image possess attributes 
such as color, or some other representation of physical data. This problem has applications 
in medical imaging where clinicians seek the best available tools to perform surgery in a 
non-invasive manner. The metamorphosis approach to this problem was formulated in Holm 
et al. [2008]. This approach can be cast into the bundle picture as follows. 

Let A/" be a manifold of deformable objects (i.e. possible images) on a manifold Q. For 
example, M = Emb (M, Q), the embeddings of a manifold M, into Q, oi M = Imm (M, Q), 
the immersions of M into Q. Suppose that the diffeomorphism group Q := DiS{Q) of Q 
acts on A/" and consider the trivial fiber bundle P = M x (A/" x Q) — M on which the 
diffeomorphism group acts on the right by the action 

{t,ri,g)h := {t^h^^ 077, /i or]). 

The projection is given by 

7rs,P : M X (TV X ^) ^ S := M X TV, 7r2,p(t, r/, g) = {t,go rj) =: {t, n). 

In the context of metamorphosis, r] is called the template, g the deformation and n = g o rj 
the image. Note that X = M and this is an example of classical Lagrangian reduction (see 
particular case iv above), for a given Lagrangian L : M x TTV x TQ —>■ M. In the context of 
metamorphosis, the Lagrangian does not depend on time and is given by 

Lig,g,v,v) = ^UWl + ^\\T9°v\\lf (5-12) 

where r G M is a parameter, ||-||g denotes a ^-invariant metric on TQ, and ||-||a^ denotes a 
metric on TTV. 

The Lagrangian in this context is interpreted as the cost of using the controls and the 
aim is to minimize 

S = L{r],r],g,g) dt, 
Jo 

where the initial image uq = g^o rjQ and the final image fii = gi o rji are given. 
The convective velocity 

A {vt, v^, Vg) = g~^Vg (5.13) 
provides a suitable connection form. Applying (2.9) yields 

Pa [it,Vt,gt,Vt,9t)]g = (Tnx,p{t'Vt,9t,Vt,9t), Kt,r]t, gt), A{t,T]t, gt,r]t, gt)}^ 
= (t,nt,ht,ut) 

where Ut = gt° gt^, rit = gt o r]t, and g = 3C{Q) is the Lie algebra of the diffeomorphism 
group Q. Note that the time derivative of nt is given by the formula fit = Ut o nt + Tgt o fji. 
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Denoting the reduced Lagrangian associated to L by / = l(t, n,h,u) : M x TAf x g ^ M the 
Lagrange-Poincare equations in this case read 



61 . 61 

= U, 

(5.14) 



dt-^ \- £ut T — 0) 
dut dut 



61 _Q 
Dt 6ht 6nt ' 

where V is a torsion free covariant derivative on JV and the fact that the group action is a 
right action is carefully noted. The reduced Lagrangian associated with L in (5.12) reads 

l{n,n,u) = ]r\\u\\l + r^\\n -uo nW^. 

Equations (5.14) are equivalent to those in Holm et al. [2008]; but they are simpler in form 
and expressed in different variables. In these variables, equations (5.14) have split into 
horizontal and vertical parts with respect to the flat connection (5.13), thereby resulting in 
their zero right-hand sides. 



5.3 Covariant formulation of metamorphosis image dynamics 

Metamorphosis image dynamics for immersions may be placed into a covariant setting. This 
is achieved by replacing X = ]RbyX = ]RxM and Imm(M, Q) by Q, and by considering 
the trivial fiber bundle 

T^x,p : P = X X {Q X g) ^ X, (t,m, q, g) t-^ (t, m) =: x. 

Now, let the diffeomorphism group Q = DiS{Q) act on P by the right action 

{x,q,g)h = {x,h~^{q),go h), 

and obtain the principal bundle 

7rs,p : P = X X [Q X g) ^ H = X X Q, {x,q,g) {x, g{q)) =: {x, n). 

In this framework, templates and deformations are sections of ttx,p where the deformation 
is further specified to be independent of the variable m G M. Concretely, for o" G F (vr^.p), 

a {t, m) = (t, m, r] (t, m) , g^t,m)) , 

where r] : X ^ Q and g : X ^ G = DiS{Q). 

The restrictions required to mimic the classical metamorphosis are: 

rjt := T] {t, ■) G Imm(M, Q) and g{t,mi) = g{t,m2), for all G M. (5.15) 

For the first restriction no constraint is necessary since our the requirement is that 

rank {T^r]t) = dimM, for all t G M, (5.16) 
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where T^rji : TM — > TQ denotes the tangent map of rjt, the variable t being considered as 
a parameter. 

The first condition in (5.15) formally defines an open subset of the space of curves with 
values in the manifold JF(M, Q) of smooth maps from M into Q. To see this, recall that 
the space Imm(M, Q) of immersions is an open subset of JF(M, Q) and that F (vrx,i;) may be 
identified with the space of curves in J-'{M, Q). 

The second condition in (5.15) is enforced by a Lagrange multiplier. Consequently, the 
new Lagrangian reads: 

L(A,A)=L(jV)-((A,(T^^^)o^-i)), 
where the Lagrange multiplier is a section 

X eT {']Tx,L{T*M,s*)) , where g* = Q^{Q), 

and ((■ , •)) denotes the pairing between the spaces L(T*M, g*) and L(TM,q). Note that 
here T^g denotes differentiation of the section dependence on M, not the argument of the 
diffeomorphism in Q. 

One creates analogous Lagrangians as before, where now the spatial derivatives of rj are 
considered to be independent variables in the theory. That is, 

J V = (x, rj, g,f,dt + T^'t], gdt + T^' g) . 

An example of a possible Lagrangian for metamorphosis is 

L{g,g,T''g,V,V,T''v) = l\m + ^2\\^^^^^ 

where || ■ ||q is norm associated to a metric on Q and || • ||q is associated to a vector bundle 
metric on L{TM,TQ). Note that in the previous formula the Lagrangian is interpreted as 
being defined on an arbitrary element of the first jet bundle J^P and not necessarily on the 
first jet extension of a section. Therefore, fj denotes an arbitrary element in TjjQ and T^^r] 
is an arbitrary element in L{TmM,TriQ). 

A diffeomorphism h E Q acts on the first jet extension j'^o", as follows: 

(x, r/, g,r]dt + T^^rj, g dt + T^' g) ^ [x, h'^ ori,goh, Th'^ o{r]dt + T^^rj), {g dt + g) o h) ; 

thus the Lagrangian L is ^-invariant. 

Fixing the connection A{VxiVq,Vg) = g^^Vg and writing {t,m) = x yields 



a{x) = [fT(x),aM]^(Q) = [(x,r7(x),c/^),c/^^(^:,. + r*VOjx(Q) 
= {x,n{x),g^g-^dt+ {T^'^g^) o g'^) , 

where n{x) = g.j.{ri{x)) which is the same definition of n as in §5.2. The Lagrange-Poincare 
equations can now be written in the convective representation from §4.1 as 

( SI ..r r 51 A n S^ 



42 



6j^n J 6n 
I {T''g) og-'=0 



Seeking a solution with 
yields the equations 



51 



6{T^g)og 



-1 



dt— + = 0, 

OUt OUt 

f 61 \ 61 ^^-^^^ 
\6j^n J 6n 

Note that the first equation above is identical to the first equation in (5.14) whilst the second 
equation (5.17) is a covariant analogue of the second equation in (5.14). 

This covariant formulation takes the problem from classical Lagrangian reduction (as 
in particular case iv), to a covariant problem that does not fall into any of the particular 
examples. 

One of the potential advantages of such a transformation would be to apply multisym- 
plectic integrators. The infinite-dimensional approach introduces special dependence on the 
time variable and is formulated on an infinite-dimensional manifold. The covariant formula- 
tion, on the other hand, takes advantage of exchange symmetry t ^ rrij for any component 
rrij of m to formulate the problem in a multisymplectic fashion. If there is no coupling in the 
Lagrangian between the fibers of AdP and J^S, then the multisymplectic problem is posed 
on a finite dimensional manifold. 



6 Conclusion and future directions 

This paper has presented a framework for Lagrange-Poincare reduction that unifies the 
approaches taken in the particular cases i - iv and extends the Lagrange-Poincare theory 
beyond the scope of those cases. On one hand, the work of Castrillon-Lopez et al. [2000] 
and Castrillon-Lopcz & Ratiu [2003] has been extended to apply to the general fiber bundle 
case. On the other hand, the classical Lagrange-Poincare theory developed in Ccndra ct al. 
[1998] and Ccndra ct al. [2001] has been extended to the field theoretic setting. 

Two surprising results have appeared that concern the integrability conditions associated 
with both reconstruction and the convective representation of the dynamics. First, the re- 
quirement of an additional condition for reconstruction first observed in Castrillon-Lopcz 
Ratiu [2003] was found also to occur here, even though less geometric structure is present. 
Second, the convective representation may not always exist for an arbitrary problem. This 
observation highlights the importance of the geometric tools used in formulating the frame- 
work. Even though the convective representation may not exist, the general AdP- valued 
objects do exist and they can be used either to study the dynamics or to find an alternative 
representation. 

Also note the large range in applications: the examples given here range from the classical 
Skyrme model in physics to molecular strand dynamics in biology, from metamorphosis image 
dynamics in computer science and control theory to the isoperimetric problem in mathemat- 
ics. This range of examples is compelling and an accurate reflection of the unifying power 
of the framework developed. 
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The covariant metamorphosis example in §5.3 has revealed an interesting property of the 
Lagrange-Poincare theory. Namely, the covariant expression of existing classical Lagrange- 
Poincare problems in the present framework produces a geometric reformulation of the prob- 
lem. Further investigation into this process of geometric reformulation could lead, for exam- 
ple, to new applications of multisymplectic integrators. 

Finally, the Kelvin-Noether Theorem has been extended to the Lagrange-Poincare field 
setting. This extension is two-fold. Firstly, the Kelvin-Noether theorem is usually stated 
for Euler-Poincare systems where there is no shape space, the Theorem now to applies to 
Lagrange-Poincare systems where a shape space is present. Secondly, the Kelvin-Noether 
Theorem now extends from the classical context to the covariant context. This result is 
particularly important, since it constitutes a major tool for gaining qualitative information 
about any problem formulated within the scope of the Lagrange-Poincare field framework. 
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